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The bosonic strictly isospectral problem for Demkov-Ostrovsky (DO) effective potentials 
in the radially nodeless sector is first solved in the supersymmetric Darboux-Witten 
(DW) half line (or /-changing) procedure. As an application, for the k = 1 class, if 
one goes back to optics examples, it might be possible to think of a one-parameter 
family of Maxwell lenses having the same optical scattering properties in the nodeless 
radial sector. Although the relative changes in the index of refraction that one may 
introduce in this way are at the level of several percents, at most, for all DO orbital 
quantum numbers I > 0, the index profiles are different from the original Maxwell one, 
possessing an inflection point within the lens. I pass then to the DW full line (or jV- 
changing) procedure, obtaining the corresponding Morse-type problem for which the 
supersymmetric results are well established, and finally come back to the half line with 
well-defined results. 
PACS number(s): 03.65.-w, 11.30.Pb 



1. INTRODUCTION 

Recently, the Schrodinger equation at zero energy with the Demkov-Ostrovsky 
(hereafter DO) class of focusing potentials 1 

V (r) = ^ (1) 

{r/Rf[{r/R)- K + {r/RYY 

has been studied in the supersymmetric quantum- mechanical context 2 . In Eq. (1) 
w > is a coupling constant, £q is an energy scale, R > is the spatial scale 
of the potential, k is the Lenz-Demkov-Ostrovsky parameter (n — 1/2 class was 
related to the atomic aufbau 3 , while the k = 1 case is a wave approach to the 
Maxwell fish-eye lens 1 ). It has been shown that the radially nodeless sector of 
the DO wave problem may be worked out in the supersymmetric approach. As 
is quite well known, the essence of Witten's supersymmetric quantum mechanics 
4 (Darboux procedure, in mathematics literature) is a pair of Riccati equations, 
which, when the particular solution is used, entails the boson-fermion symmetry 
at the level of standard one-dimensional Schrodinger equations. Moreover, when 
use is made of the general Riccati solution one can solve the isospectral problem, 
i.e., one is able to generate, in the most simple procedure, classes of one-parameter 
families of potentials having the same eigenvalues as the original potential. The 
one-parameter family of fermionic DO potentials for the zero energy case has been 
already obtained in 2 . 
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The organization of the paper is as follows. In Sec. II, I shall obtain the "half 
line" one-parameter family of bosonic DO effective potentials at zero energy, in the 
nodeless sector, having the same fermionic superpartner. Moreover, I shall focus on 
the optical application, introducing, in the nodeless radial sector, a one-parameter 
family of generalized Maxwell indices of refraction, in the sense that they are defined 
for each partial wave and depend on the family parameter A. Then, in Sec. Ill, I 
shall use the Langer transform to go to the full line Schrodinger equation for which 
the Darboux-Witten (DW) procedure applies in a standard way and return to the 
half line with simple yet well-defined results. A conclusion section ends up the 
paper. 

2. DO-ISOSPECTRAL PROBLEM ON THE HALF LINE 

I shall use the double DW procedure to first delete and then reintroduce a 
normalizable state in a spectrum, taking the initial DO Schrodinger equation at zero 
energy, which for the quantized values wn,k — (2k) 2 [N + (2k)~ 1 ][N + (2k) _1 — 1] 
is known to have a zero-energy, regular, normalizable state whose radial part is 
(P = r/R) 

R <n\ - ^ Nl()l r (2l+l)/2K+l/2 (c) 

Hm yP> - ( 1 + /0 2 K )( 2 ;+i)/2 K G JV-i-i/« Is) ' W 

where N = n+ — 1)/, n = n r + 1 + 1, n r = 0, 1, 2, arc the "total", "principal" 
and "radial" DO quantum numbers, respectively, I is the DO orbital one, C£(£(p)) 
are Gegenbauer polynomials of variable £ = (1 — p 2ft )/(l + p 2li ), and A/jvz are the 
normalization constants that I shall not take into account henceforth. The degree 
of degeneracy of the zero-bound state is N 2 similar to the hydrogen atom. In the 
n r = sector the Gegenbauer polynomials are unity for any degree p, so they do 
not enter the following. 

To generate the bosonic isospectral family I need the general solution of the 
fermionic Riccati equation 

^/ = +^ + w2M. (3) 

for which the particular solution is the following 2 

W «(P) = " - nlL > ( 4 ) 

p p(1 + p 2k ) 

and where is the fermionic effective superpartner 2 . Substitution of the general 
solution supposed to be of the form W K = V~ 1 + W K in Eq. (3) leads to the following 
simple differential equation for V K 

-^ + 2W K V K = -l (5) 

with the general solution (the particular superpotential is the negative logarithmic 
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derivative of the radial factor f K ) 

K,a = /- 2 (a + £ f K {p)dp^j , (6) 

where A is a positive, real parameter, and the DO radial factor in the n r = sector 
is 2 

p i+i 

f K (ti = (1 +p 2 K )(2;+l)/2 K • ( 7 ) 

Denoting the integral in Eq. (6) by Io- K (p) the general superpotential reads 

W A , K = W K + ^ln[/ 0;K (p) + A] , (8) 

as one can show easily. Then, the isospectral effective potential family can be 
written in terms of the original effective DO potential and the radial factor f K as 
follows 

U bos , K (p; A) = U-(p) 2[ln(7 0;K + A)]" = U~(p) + ^ , (9) 

the primes denoting derivatives with respect to p. It remains to calculate the 
integral Iq- k in order to consider the problem as solved. For that, the trigonometric 
representation is the most convenient. With p = [tan/?] 1 /* the integral turns into 



Io; K = - f [cos/3] "[sin/3] + « 
K Jo 



d/3 . (10) 



Then, for the physical 1/2 and k = 1, one can use the formulas 2.524.1 

and 2.518.1, respectively, in 5 , as the most compact ones, in order to get 



, [8m0\«+* ™ k+1 (2l + 2)\ [cos/?p-3 

7o - = 2 i i^r d ^ 2 t { - l) m + 2-ky. [2* -3] UL) 

and 

/ . , _ 

/ [cos/3] 2 cos/3 



,,.,= , I^^.BEfl^i.pB + ^-afl-*. (12) 



Finally, for the last integral in Eq. (12) I used the formula 2.513.1 in 5 leading to 

_ [sin/3] 2 '+ 1 (2/ + l) (2») (-1)'^ (2Q! [sin(2Z - 2fc)/3] 

0;1 " cos/3 4' G!) 2 ^ ^ <+1J 4' ^ j fc!(2Z-fc)! [2Z - 2k] 

(13) 

2.1. Application 

The interesting physical point may be seen more directly in the Maxwell fish- 
eye case, k = 1. Notice that now N = n. Suppose we use an optical physics 
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example. Then, one might think of a whole family of indices of refraction of the 
Maxwell type, i.e., a family of lenses with the same remarkable optical properties. 
Since after performing the double DW procedure the centrifugal barrier is the same 
for the whole bosonic family and equal to that of U~^, 6 following Demkov and 
Ostrovsky 1 , to the bosonic family of Maxwell fish-eye potentials 

(2Z + lj(2Z + 3) 4/i/; 2ft 

Vlip ' l ' A) - (1 + p 2 ) 2 " J^Tx + Jf^+W (14) 

would correspond the family of Maxwell indices 



m{ P ; I, A) cx y/-Vi(p;l,X) . (15) 

The one-parameter family of Maxwell lenses that might be thought of in this way 
are identical in the sense that, quantum-mechanically, the scattering matrices are 
identical 6 . I emphasize that all these n(p; I, A) are valid in the n = I + 1 sector. (See 
Figs. 1 and 2.) 

It is known in the literature 6 that pair combinations (including double-ones, 
as the double DW used here) of all the three known procedures to delete and 
reintroduce a nodeless, normalizable state in a spectrum, i.e., the Darboux-Witten, 
the Abraham-Moses, and the Pursey ones, generate five distinct, one-parameter, 
"isospcctral" families of (effective) potentials. However, the supplementary four are 
not exactly isospectral. Nevertheless, their quantum-mechanical scattering matrices 
are related by simple relationships 6 . 

Since I have used the double DW procedure I first deleted and then reinserted 
the zero-energy bound state, or a sector of it. The price for such an action is the 
damping of the new radial factors as follows 

hosAp) = ■ (i6) 

From the plots I did, one can notice that the "physical" effect of the family param- 
eter A on the index of refraction, though not important - it is at the level of several 
percents and goes down strongly with increasing I - is sufficient to generate an in- 
flection point in the index profile. Its physical nature should be further examined. 
On the other hand, A is acting as a huge damping parameter for the radial factor. 

3. FROM THE HALF LINE TO THE FULL LINE AND BACK 

What I have discussed so far has to do with the half-line (radial) problem. It 
is a supersymmetric construction which for the very close hydrogen problem has 
been questioned by Haymaker and Rau 7 . This is so because the interpretation of 
the DW construction along the lines of Fadeev 8 , Sukumar 9 , or Kostelecky and 
Nieto 10 , i.e., as changing the angular momentum, is a delicate point, and actually 
an important open issue in Physics. Indeed, to associate the DW procedure with a 
change n — > n — 1 at fixed I as suggested by Haymaker and Rau seems more natural. 
Therefore, I shall tackle this point as well. 
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The idea is to pass from the half line case to the full line one by the Langcr change 
of variable x = hip and wave function 4>{x) = cxp(— x/2)u. Then the standard DO 
Schrodingcr equation at zero energy turns into the Rosen-Morse (RM) problem 

d 2 , , 1,2 (n-|)(n+i) 
dx* +[n 2 ] cosh 2 * 

where I made use of I + 1 = n, considered as the lowest eigenvalue and fixing the 
zero of the energy scale on it. I shall also make the shift n& = n + 1/2, taking only 
the integer part [rib] , since I am not interested in the WKB analysis and I need the 
standard form of the potential 

K- = _!^+l). (18) 

cosh x 

Then, Eq. (17) is just the RM eigenvalue problem with [rib] bound states, of eigenval- 
ues — [rib] 2 . The particular superpotential is W — rib tanhx and the superpartner is 
V+ = -S±iS±fll and therefore has one bound state less than V . The supersymme- 
try of the full line RM potential is well established, including the strictly isospectral 
issue, see, e.g., 11 . Multiple-parameter families of strictly isospectral — cosh -2 x po- 
tentials can be obtained easily and a nice picture of the translational motion of the 
wells, each with one bound state, towards ±oo as the family parameters are going 
to the Abraham-Moses limit (— 1~) and Pursey limit (0 + ), respectively, is also well 
known. For example, in the case of a single bound state, the one-parameter family 
of — cosh -2 x potentials reads 

Vo(X ) = -2cosh- 2 ( a; + iln(l + l/A )) (19) 

where one can see the parameter imbedded in the argument of a logarithm, adding 
to the argument of the hyperbolic function, which is a general feature of this case 
and the reason for the translation of the potentials. 

Exactly as suggested by Haymaker and Rau, after all the supersymmetric job 
on the full line, one should go back to the half line to obtain the right superpartner 
and isospectral families. By doing so, one gets the DO particular superpotential 
W = (| — n)^(p), a superpartner effective potential U + = 1(1 + l)/p 2 - (21 + l)(2l — 
1)/(1 + p 2 ) 2 , whereas the effect of the isospectral parameter shows up in a rescaling 
of the p variable of the type p\ = \Jl + j^Poc i where is the initial coordinate, 
corresponding to Ao = oo. In other words, we have a rescaling of the radius of 
the fish-eye lens occuring simultaneously with the change n — ► n — 1. The fish-eye 
lens of radius R\ = Ry / X /(\n + 1) has a zero-bound state of degeneracy reduced 
by one unit with respect to the lens of radius R. This is a counterpart of what 
happens in the hydrogen case, where one has the simultaneous change n — > n — 1 
and Z^Z(1- 1/n). 7 

In the atomic aufbau case the same results hold with only minor changes in the 
RM problem, which this time turns to apply to the potential — iV(iV+l)/4cosh (x/2) 
of eigenvalues — N 2 /4, where N = 21 + 1. 



(17) 
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4. CONCLUSION 

In conclusion, I presented a rather detailed discussion of the DO wave problem 
in the nodeless radial sector by employing DW tools, with special focus on the fish- 
eye lens case. While the ^-changing half line procedure does not look very natural, 
the iV-changing full line one appears to be more acceptable. Of course, it is more 
convenient not to face the centrifugal singularity by sending it to — oo, avoiding all 
related difficulties. One just pass in this way to an easy to handle rcficctionless and 
shape invariant RM problem. The centrifugal singularity is not "active" . Moreover, 
the isospectral parameter plays different roles, as a strong damping in the half line 
case, and as a "dilatation" parameter when the full line procedure is used. 

As an important by-product of this study I introduced generalized 1- and A- 
dependent families of Maxwell indices of refraction, which presumably might be 
related to surface-scattering optical phenomena. Indeed, the squares of the radial 
factors are peaked close to the surface p = 1. Finally, as one might guess from the 
above analysis, there might be potential applications of DW techniques in graded 
index optics. 
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Figure Captions 

Fig. 1. The original Maxwell refractive index um, the "isospectral" one rii SOl 
their relative ratio — 1, and the square of the radial factor / [Eq. (16)] for (a) 
I = 1, A = 1 and (b) I = 1, A = 10, respectively. 

Fig. 2. The same functions for (a) I = 2, A = 1 and (b) I = 2, A = 10, 
respectively. 

In both figures, I have used n iso /riM = 1 + where V\ is the negative 

of the A-depending part in Eq. (14), whereas Vm is the negative of the first term 
therein. Since V\ may be negative, the relative ratio may be negative as well. The 
indices have been normalized to (I + 1/2) 2 to be as close as possible to the Maxwell 
normalization constant n(0) = 2, and also to include the case I = 0. 
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